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ABSTRACT 

In  an  earlier  thesis,  Ellsworth  showed  that  the  radiated  power  output  (quantified  by 
the  gain  in  radiation  resistance  presented  to  the  transducer)  and  directivity  of  a  compact 
underwater  transducer  can  be  significantly  improved  by  the  use  of  a  system  of  resonant 
scatterers.  These  resonant  scatterers  were  termed  "sympathetic  resonators."  In  the 
present  work,  we  verify  the  previous  findings  and  extend  the  work  to  additional  cases  of 
interest.  Specifically,  we  (1)  employ  equations  requiring  fewer  assumptions  and 
implement  them  using  more  accurate  numerical  techniques,  (2)  reproduce  Ellsworth's 
calculations  and  provide  graphical  results  for  conical  and  planar  circular  configurations 
for  the  case  of  six  passive  acoustic  resonators,  (3)  extend  the  theory  to  and  provide 
results  for  two  typei  of  linear  arrays  of  passive  acoustic  resonators,  (4)  and  finally,  we 
consider  an  acoustic  Yagi-Uda  antenna-like  configuration  of  passive  acoustic  resonators 
to  improve  the  directivity  of  an  underwater  transducer.  The  results  indicate  that  the 
conical  configuration  of  six  resonators  is  the  best  choice  in  terms  of  the  gain  in  radiated 
power  output,  directivity  along  an  acoustic  axis,  and  minimum  required  number  of 
resonators. 
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I.    INTRODUCTION 

A.      BACKGROUND 

It  has  been  shown  that  a  coupled  set  of  resonant  scatterers  will,  at  frequencies  near 
the  resonance  of  a  single  scatterer,  exhibit  a  quasiresonant  condition,  resulting  in  a 
significantly  higher  scattering  cross  section  than  a  single  scatterer  [Ref.  1],  [Ref.  2]. 
This  phenomena  has  been  shown  by  I.  Tolstoy  to  provide  significant  gain  in  the  scattered 
acoustic  pressure  to  a  plane  wave  when  incident  upon  a  system  of  scatterers  as  in  Figure 
1-1  [Ref.  3].  Tolstoy  demonstrated  that,  particularly  for  the  line  array,  the  amplification 


Plane  Wave 


Resonant  Scatterers 

Figure  1-1    Plane  wave  incident  upon  a  system  of  resonant  scatterers. 

of  the  scattered  wave  is  sensitive  to  the  direction  of  the  incoming  plane  wave  [Ref.  3]. 
In  his  thesis  [Ref.  4],  J.  Ellsworth  considered  the  effect  of  the  presence  of  a  system  of 
resonant  acoustic  scatterers  on  a  transducer  in  the  coupled  acoustic  field  as  in  Figure  1-2. 


Ellsworth  referred  to  these  resonant  scatterers  as  "sympathetic  resonators."  We  will 
refer  to  these  scatterers  as  passive  acoustic  resonators.  Two  geometries  were  previously 
considered.  The  first,  called  the  planar  configuration,  placed  the  transducer  in  the  center 
of  a  circular  array  of  resonators.    The  second,  called  the  conical  configuration,  is  the 


Resonant  Scatterers 


Transducer 

Figure  1-2   A  transducer  in  a  coupled  field  of  resonators. 

same  as  the  planar  configuration  except  that  the  transducer  is  displaced  one-quarter 

1 

wavelength  out  of  the  plane  of  the  resonators.  The  results  demonstrated  a  gain  in 
radiation  resistance  presented  to  the  transducer.  This  means  that,  for  a  displacement 
limited  transducer,  more  power  can  be  radiated  into  the  medium.  Ellsworth  also 
provided  beam  patterns  demonstrating  that  the  transducer  will,  in  the  presence  of  passive 
acoustic  resonators,  exhibit  directionality  along  an  acoustic  axis. 

B.      OBJECTIVES 

The  objectives  of  this  thesis  are  to  (1)  refine  the  theory  presented  by  Ellsworth 
using  equations  requiring  fewer  assumptions  and  implementing  those  equations  with  more 
accurate  numerical  methods,  (2)  reproduce  Ellworth's  work  using  the  more  accurate 


techniques  and  provide  three-dimensional  plots  for  the  case  of  six  resonators  (the  six 
resonator  results  were  requested  by  the  Naval  Ocean  Systems  Center),  (3)  develop  three 
dimensional  beam  patterns  for  all  cases  of  interest,  (4)  develop  the  theory  for  and  provide 
results  for  two  types  of  line  arrays  of  passive  acoustic  resonators,  and  (5)  to  consider  an 
acoustic  Yagi-Uda  [Ref.  5]  antenna-like  configuration  of  scatterers  to  improve  the 
directivity  of  an  underwater  transducer. 

Our  measures  of  effectiveness  in  the  above  work  will  be  (1)  the  gain  in  radiation 
resistance  presented  to  the  transducer  with  the  system  of  passive  acoustic  radiators 
present,  (2)  three-dimensional  plots  of  that  gain  versus  the  system  parameters,  (3)  the 
directivity  index  of  the  transducer  in  the  presence  of  passive  acoustic  radiators,  and  (4) 
the  shape  of  and  characteristics  of  the  far-field  beam  pattern  for  each  configuration 
considered.  In  support  of  these  objectives,  programs  were  written  in  Mathematica  [Ref. 
6]  (a  system  for  doing  mathematics  oy  computer)  on  Macintosh,  IBM,  and  SUN 
workstation  computers.    An  example  program  is  provided  in  the  appendix. 

The  remainder  of  this  thesis  is  organized  as  follows.  We  first  develop  the  theory 
for  a  transducer  in  a  coupled  acoustic  field  with  passive  acoustic  resonators.  We  provide 
results  for  the  two  circular  configurations,  planar  and  conical.  We  analyze  and  provide 
results  for  two  types  of  line  arrays  of  resonators,  termed  line  type  I  and  line  type  II.  We 
consider  the  acoustic  Yagi-Uda  array  of  resonators.  Finally,  we  discuss  the  results  and 
provide  conclusions  and  recommendations. 


D.  THEORY 

A.      THE  NETWORK  MODEL 

We  will  develop  the  network  model  for  a  coupled  network  of  a  transducer  and  N 
passive  acoustic  resonators  by  considering  each  element  as  a  lumped  acoustic  system. 
The  development  is  similar  to  that  presented  by  Ellsworth  [Ref.  4].  We  will  represent 
each  lumped  acoustic  element  as  an  acoustic  impedance,  choosing  pressure  and  volume 
velocity  as  mechanical  variables.  In  Section  1  we  develop  the  network  representation  of 
a  single  resonator  in  an  acoustic  field.  In  Section  2  we  develop  the  network 
representation  of  a  transducer  in  an  acoustic  field.  In  Section  3  we  combine  the  results 
of  Sections  1  and  2,  and  develop  a  network  model  for  a  transducer  coupled  to  N 
resonators  by  an  acoustic  field.  We  allow  for  interacions  between  each  element  and 
every  other  element.  This  network  will  be  used  to  analyze  systems  of  passive  acoustic 
resonators  in  the  chapters  that  follow.  Finally,  in  Section  4,  we  discuss  the  required 
input  parameters  for  the  model  developed. 

1.       A  RESONATOR  IN  AN  ACOUSTIC  FIELD 

Figure  2-1  shows  a  lumped  acoustic  model  of  a  single  resonator  in  an  acoustic 
field.  The  network  equations  for  this  system  are: 

P,    -    -Z'U,  C-la) 

and 


n     Ux     +    Px 
where  px     =       acoustic  pressure  at  the  surface  of  the  resonator 


(Mb) 


Resonator 


Acoustic  Field 


Figure  2-1    Lumped  acoustic  model  of  a  resonator  in  an  acoustic  field. 


U]    =       acoustic  volume  velocity  of  the  resonator 
Zf    =     acoustic  impedance  of  the  resonator 

Z°[    =     acoustic  radiation  impedance  presented  to  the  resonator  by  the 

medium  with  all  other  active  or  passive  elements  blocked 
pjb    =      the  acoustic  pressure  acting  on  the  surface  of  the  resonator  when 

it  is  blocked  (i.e.  JJX  =  0) 
D     =      diffraction  constant  [Ref.  7] 

pf    =      free-field  pressure  incident  on  the  resonator,  i.e.  the  pressure  that 
would  exist  at  the  location  of  the  resonator  if  it  were  not  present. 


The  free-field  pressure  can  be  found  from 

pf    =    £  Zg  U  (2-2) 

F\  m         \m         m 

where  m  increments  over  all  other  active  and  passive  elements  present  and 

ZU?    =     ^e  ft®e-field  acoustic  transfer  impedance  from  element  m  to  the 

resonator 
Um    =      acoustic  volume  velocity  of  element  m. 

Bobber  [Ref.  7]  defines  the  diffraction  constant  as  the  following  ratio 


D,    -    El.      .  (2-3) 

pf 


The  blocked  pressure  is  the  average  pressure  over  the  surface  of  the  resonator  if  its 
volume  velocity  is  zero.  The  surface  pressure  is,  in  general,  the  sum  of  the  blocked 
pressure  and  the  radiated  pressure,  as  indicated  in  Figure  2-1.  The  radiated  pressure  is 
given  by 

/>,"*    =    Ux  Z°[  (2-4) 

and  so  the  total  pressure  at  the  surface  is 

p,   -  ri  +P?  -  dp*  *  v&; .  (2-5) 

For  notational  convenience,  we  will  include  the  diffraction  constant  as  part 


of  the  acoustic  transfer  impedance  between  two  elements.     Let  zar  be  the  transfer 


impedance  from  element  m  to  element  n.    We  can  then  write 


.  E 


m 


Z     U 

**nm    um 


(2-6) 


The  radiation  impedance  is  a  function  of  the  properties  of  the  medium  and 
the  physical  size  and  shape  of  the  object  radiating  into  the  medium.  The  input  acoustic 
impedance  is  a  function  of  the  properties  of  the  resonator. 

2.       A  TRANSDUCER  IN  AN  ACOUSTIC  FIELD 

Figure  2-2  is  the  lumped  acoustic  model  of  a  transducer  in  an  acoustic  field. 
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Figure  2-2   Lumped  acoustic  model  of  a  transducer  in  an  acoustic  field. 
The  network  equations  for  this  system  are 


T  I0  ~  Z"  U0 


(2-7a) 


Zar    T ,  b 

oo  U0  +  p0 


(2-7b) 


and 


V0     =    ZEB  I0-TU0  (2-7c) 

where         p0    =       acoustic  pressure  at  the  surface  of  the  transducer 

U0    =       acoustic  volume  velocity  of  the  transducer 

Z20    =      open-circuit  acoustic  impedance  of  the  transducer 

T       =     transduction  coefficient  of  the  transducer 

Zar 

00    =     acoustic  radiation  impedance  presented  to  the  resonator  by  the 

medium  with  all  other  active  or  passive  elements  blocked 
p*    =      the  blocked  acoustic  pressure  at  the  surface  of  the  transducer 

ZEB    =      blocked  electrical  impedance  of  transducer 
V0     =      voltage  across  the  electrical  terminals  of  the  transducer 
10     =       current  through  the  electrical  terminals  of  the  transducer. 
The  blocked  acoustic  pressure  is  found  as  before.    The  electrical  quantities 
and  physical  characteristics  of  the  transducer  will  not  be  considered  further  here  but  are 
included  for  completeness.   The  transducer  is  taken  to  be  compact  and  we  are  interested 
in  the  radiation  impedance  presented  to  it.    The  volume  velocity  of  the  transducer  will 
be  held  constant  (assumed  to  be  displacement-limited).    For  a  displacement-limited 
transducer,  the  power  radiated  can  only  be  increased  by  increasing  the  radiation  loading. 


3.       A  TRANSDUCER  IN  THE  PRESENCE  OF  PASSIVE  ELEMENTS 

Figure  2-3  is  the  network  model  of  N  interacting  passive  acoustic  resonators 
coupled  to  a  single  transducer  by  an  acoustic  field.   This  model  combines  the  results  of 
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Figure  2-3  Network  model  of  N  interacting  passive  acoustic  resonators  and  a  transducer 

the  previous  two  sections  and  allows  for  interactions  between  each  element  and  every 
other  element.   The  network  equations  are 


p0    =    TI  -  Z"  U0 


(2-8a) 


and 


p„    =    -z;  Un  n-0,1,2, N  <2"8b) 


p.    -    £  Z-  P.  B-1,2,3 W  (2-8c) 

where  #    =        number  of  passive  acoustic  resonators 

pn    =       acoustic  pressure  at  the  surface  of  element  n 

Un    =      volume  velocity  of  element  n 

Z°rn    =      acoustic  radiation  impedance  seen  by  element  n  with  all  other  active 

or  passive  elements  blocked 
Zar    =     acoustic  transfer  impedance  from  element  m  to  element  n. 

We  are  interested  in  the  total  acoustic  radiation  impedance  seen  by  the 
transducer  in  the  presence  of  the  passive  acoustic  resonators.  This  quantity  can  be 
defined  as  follows 

Z,l,     -    £.  (2-9) 

Equation  2-8  can  be  written  for  the  case  of  N  =  6  resonators  as  the  following 
matrix  equation 
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In  order  to  simplify  the  notation,  we  will  write  equation  2-10  as 


/'o 


'00 


{zs} 


n 


{zs} 


E/„ 


W 


(2-11) 


where 


W 

M 


is  a  lxN  column  matrix 


is  a  NxN  square  matrix 


and  the  superscript  "T"  indicates  the  transpose  of  a  matrix. 


In  order  to  find  z,^;  , first  we  replace  fo\  by\-2ajj  }  and  normalize  by  U0, 


resulting  in  the  following  equation: 


yar 
^loial 

Zoo               1  m0J 

». 

_{zz)  [z;  ♦  z;  8^ 

u, 


rel 


(2-12) 
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where  5       =    l\         n  "  m 

B"1  0         n  itm 


and  Urme!    =    Hk 


We  can  represent  the  matrix  of  equation  2-12  as  the  following  system  of  equations: 

zz*  -  zz  *  WKi  (2_13) 


and 


{0}     =    {Z-}     ♦    [ZZ  +  ZBfl  5 J  {u?}  n,m  =  \,2,...,N  .         (2~14) 

Solving  for  I^C  j  in  equation  2-14  yields 

{ur:!}  -  -  [z-  +  z;  aj"1  {zj} .  (2"15) 

Equation  2-15  gives  the  relative  volume  velocities  of  each  resonator  with  respect  to  the 
transducers  volume  velocity.    We  substitute  equation  2-15  into  equation  2-13  to  obtain 

the  following  expression  for  Z^,: 


=     7l 


'tola!  *-00 


-  [z:]T  &  +zn°  5J"  {25}  •  (2_16) 


The  total  radiation  impedance  presented  to  the  transducer  will  be  used  to  find  the  gain 
in  radiation  resistance  with  resonators  present. 
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4.       THE  INPUT  PARAMETERS 

The  physical  parameters  that  we  must  provide  in  the  model  developed  above 
are  (1)  the  acoustic  impedance  of  each  resonator,  (2)  the  transfer  impedance  of  each 
element  in  the  presence  of  other  elements  and  with  the  other  elements  blocked,  (3)  the 
physical  characteristics  of  the  medium  (hydrostatic  pressure,  sound  speed,  etc.)-  The 
required  equations  are  developed  later  in  this  chapter.  The  equations  for  the  measures 
of  effectiveness  chosen  to  evaluate  the  resonator  systems  are  developed  in  the  sections 
that  follow. 

B.       GAIN  IN  RADIATION  RESISTANCE 

The  total  power  radiated  into  the  water  by  a  simple  source  is  [Ref.  8] 


n  =  y  \T° Re[f] Re[u^dt  (2_17) 


where  f    =         force  of  the  fluid  on  the  source 

u0    =       particle  velocity  at  the  surface  of  the  source 

T0    =       period  of  harmonic  motion. 

In  terms  of  the  variables  used  in  this  paper,  the  total  power  delivered  to  the  acoustic  field 
by  a  transducer  in  the  presence  of  N  passive  resonators  is  given  by 

H     =    \  U2o  Re[ZaJJ  .  (2-18) 

The  gain  in  power  delivered  to  the  acoustic  field  due  to  the  presence  of  the 
resonators  can  be  defined  as  [Ref.  4] 
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c  .  total  power  delivered  with  the  resonators  present  (2-19) 

total  power  delivered  without  the  resonators 

For  a  compact  source,  the  radiated  acoustic  power  in  the  absence  of  the  resonators  can 
be  approximated  by 

J],     -    \  Vl  Re[Z£)  .  (2-20) 

Then  we  can  write  equation  2-19  as 


Gain    =    Re[Z^  .  (2-21) 


We  will  refer  to  this  quantity  as  the  gain  in  radiation  resistance,  or  simply  gain.  The 
gain  in  radiation  resistance  is  an  important  measure  of  effectiveness  because  it  quantifies 
the  additional  radiated  power  due  to  the  presence  of  passive  acoustic  radiators.  We  will 
use  this  quantity  to  evaluate  the  benefit  of  using  passive  acoustic  radiators  in  various 
configurations. 

C.      BEAM  PATTERNS  AND  DIRECTIVITY 

The  directivity  of  a  transducer  due  to  the  presence  of  passive  acoustic  resonators 
can  be  determined  by  analyzing  the  far-field  beam  pattern.  To  find  the  beam  pattern  we 
find  the  radiated  pressure  due  to  each  element  at  far-field  points  in  various  directions. 
For  a  compact  source,  we  can  calculate  the  field  pressure  due  to  a  single  element  as 
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pT    -     Un  ZZ  QJL  iOl  (2-22) 

where  p^eU    =    acoustic  pressure  at  a  field  point  due  to  element  n 

a     =       radius  of  element  n 

n 

rn    =       distance  to  the  field  point  from  the  acoustic  center  of  element  n. 
The  total  pressure  at  the  field  point  due  to  the  transducer  and  system  of  resonators  is  then 

w=0 

where  p^M(r,6,4>)    =   the  pressure  at  the  field  point  specified  by  the  spherical 

coordinates  (r,6,<f))  due  to  element  n. 

In  the  plots  which  will  be  presented,  the  above  far- field  pressure  is  normalized  by 
the  pressure  at  the  same  point  due  to  the  transducer  alone.  What  is  plotted,  then  is  the 
following  function,  termed  BP(0,<£): 


BPiBM  -  Jf™, 


rfVA*) 


(2-24) 


We  will  plot  the  beam  pattern  defined  by  equation  2-24  in  three  dimensions  in  order  to 
easily  see  what  is  the  effect  of  the  resonators  on  the  pressure  field.  We  will  also  plot  the 
beam  pattern  on  a  log  scale  in  the  form  of  a  polar  plot.    Polar  plots  are  the  traditional 
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method  of  presenting  a  beam  pattern.  In  this  case  0  dB  will  represent  the  sound  pressure 
level  of  the  transducer  with  no  resonators  present. 

A  measure  of  how  well  the  system  projects  power  along  an  acoustic  axis  is  the 
directivity  index.   The  directivity  index  is  defined  as  [Ref.  8] 


Dl    =    10  logl0 


Plxif) 


(2-25) 


where         PaJf)    =  acoustic  pressure  due  to  the  directional  source  at  a  range  r 

along  the  acoustic  axis 
ps{r)    =  acoustic  pressure,  at  a  range  r,  due  to  a  spherical  source 

radiating  the  same  acoustic  power  as  the  directional  source. 
In  terms  of  our  equations,  the  directivity  index  can  be  show  to  be 

Dl    =    20  log10  (BPiS^J)  -  201og10(v/G^)  (2_26) 


where  Bm    and    4>m  indicate  the  spherical  angles  along  the  acoustic  axis. 

D.      GAS-FILLED  BUBBLES  AS  PASSIVE  ACOUSTIC  RADIATORS 

Tolstoy  and  Ellsworth  both  considered  air-filled  bubbles  as  passive  acoustic 
resonators  in  their  papers.  Here,  we  will  consider  air-filled  bubbles  at  a  hydrostatic 
pressure  PQ=  1  atm  (just  below  the  surface).    We  will  assume  a  sound  speed  in  water 
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of  c  =  1500  meters  per  second  and  a  density  in  water  of  p0  =  998  kg/m3.  The  ratio  of 

specific  heats  for  air  is   7  =  1.4. 

In  addition  to  the  above  constants,  we  require  the  acoustic  impedance  of  each 
resonator  and  the  transfer  impedances  between  all  elements.  The  required  impedances 
for  air- filled  bubbles  are  developed  below. 

The  mechanical  stiffness  of  an  air  bubble  in  water  is  [Ref.  8] 

s    =    12  r  a  7  P0  (2-27) 

where  a  is  the  radius  of  the  bubble.    This  stiffness  can  be  represented  as  an  acoustic 
impedance  as  follows 

Z;    -    -jllll  (2-28) 

S  c  k  a 

where  s    =        surface  area  of  the  bubble 

jfc    =         wavenumber  in  water. 
We  have  only  included  the  adiabatic  stiffness  of  the  air,  since  the  effects  of  surface 
tension,  non-adiabatic  expansion  and  contractions,  and  losses  due  to  viscosity  can  be 
neglected  for  bubbles  larger  than  1  mm  [Ref.  9]  and  frequencies  below  40  kHz  [Ref.  8]. 
The  resonant  frequency  of  a  bubble  in  water  is  [Ref.  8] 


w0     =     - 


3  7  ^0  (2-29) 
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We  will  prefer  to  consider  dimensionless  quantities  such  as  the  wavenumber  times 
distance.    Equation  2-29  can  then  be  expressed  as 


ka     = 


3  y  po  (2-30) 


PQ    C2 


For  the  values  specified  above,  the  resonant  value  of  ka  is  kar  =0.01377. 
The  acoustic  self-impedance  seen  by  element  n  is  [Ref.  8] 


P<£         ka 

22. "  eJ 


where  cot(0)     =    ka  .    This  can  be  easily  rewritten  as 


(2-31) 


P£JK_  (2.32) 

S   l+jkan 


The  transfer  impedance  from  element  m  to  element  n  is 

ka 
Zl    -    CD.^-'y  (2-33) 

where  /      is  the  distance  between  the  acoustic  centers  of  elements  m  and  n.     The 

nm 

diffraction  constant  for  a  sphere  of  radius  an  is  [Ref.  7] 


J*. 

D.    =    J. (2-34) 


\+jkan 
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We  have  taken  the  transducer  radius,  a0,  to  satisfy  ka0  =  0.2.  This  is  a  typical 
value  for  a  low  frequency  transducer.  Later  we  will  show  that  the  results  are  not 
sensitive  to  this  choice  over  a  large  range  of  values.  This  value  of  transducer  radius  was 
used  for  all  calculations  unless  otherwise  noted. 

E.      OTHER  PASSIVE  ACOUSTIC  RADIATORS 

The  network  model  developed  above  is  general  and  can  be  used  for  a  variety  of 
passive  acoustic  radiators.  Any  compact  resonator  that  can  be  characterized  by  a 
stiffness  alone,  such  as  an  air  bubble  in  water,  can  use  the  impedances  developed  for  the 
bubble  with  the  appropriate  choice  of  kar.  In  addition  to  the  air  bubble  at  1  atm  ambient 
pressure,  we  have  chosen  to  perform  calculations  for  a  bubble-like  resonator  with  a  value 
of  kar  =  0.1  (which  corresponds  to  an  air  bubble  at  50  atm  ambient  pressure).  We  can 
expect  that  a  very  thin-shelled  object  made  of  steel  or  aluminum,  an  air  balloon,  or  air 
bladder  might  be  built  with  these  characteristics.  We  have  not  attempted  to  devise  or 
construct  such  a  system,  but  rather  present  results  to  determine  whether  one  should  be 
built.  We  will  compare  the  results  for  air  bubbles  at  1  atm  ambient  pressure  and 
resonators  with  kar  =  0.1. 
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ffl.    THE  CIRCULAR  CONFIGURATIONS 

In  this  chapter,  we  present  the  results  for  the  two  circular  configurations  previously 
analyzed  by  Ellsworth  [Ref.  4].  The  tables  of  results  presented  below  were  created  using 
more  accurate  numerical  techniques  than  those  used  by  Ellsworth.  Also,  the  equations 
used  require  fewer  assumptions.  Specifically,  (1)  we  use  an  exact  equation  for  the 
radiation  impedance  from  a  sphere  (rather  than  an  approximate  formula),  (2)  we  take  into 
account  the  diffraction  constant  relating  the  blocked  pressure  to  the  incident  pressure 
wave,  (3)  we  used  a  gradient  search  method  to  find  the  maximum  gain  as  a  function  of 
the  free  parameters  rather  than  reporting  the  maximum  found  over  the  coarse  grid  of 
parameter  values,  (4)  we  calculate  the  directivity  index  for  each  case  considered,  and  (5) 
we  plot  three-dimensional  beam  patterns  in  order  to  better  visualize  the  effect  of  passive 
acoustic  resonators  on  the  acoustic  field.  Ellsworth  provided  plots  of  gain  in  radiation 
resistance  and  beam  patterns  for  the  case  N=4.  In  this  thesis  we  present  graphs  for  the 
case  N  =  6.  It  will  be  shown  that  the  six  resonator  case  is  superior  to  the  four  resonator 
case,  and  is  in  fact  the  best  choice.  In  the  following  sections  we  present  results  for  the 
planar  and  the  conical  configurations. 
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A.      THE  PLANAR  CONFIGURATION 


1.       THE  GEOMETRY 


Figure  3-1  shows  the  geometry  of  the  planar  circular  configuration  for  the 
case  of  six  passive  acoustic  resonators.  The  radius  of  the  circle  is  R  and  the  transducer 
is  in  the  plane  of  the  resonators  (x-y  plane).  We  are  interested  in  dimensionless 
distances,  so  the  radius  of  the  circle  will  be  referred  to  in  terms  of  wavenumber  times 
the  radius,  kR.  The  dimensionless  radius  of  each  resonator  is  ka  and  the  dimensionless 
radius  of  the  transducer  is  ka^  .     We  assume  a  value  for  the  transducer  radius  of 


Figure  3-1  Geometry  for  the  planar  configuration  with  six  resonators.  The  plot  is  drawn 
to  scale  for  resonators  with  ka  «  0.1. 
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ka0=0.2.  Later  we  will  demonstrate  that  the  results  are  independent  of  the  value  chosen 
over  the  range  of  interest. 

2.       GAIN  IN  RADIATION  RESISTANCE 

Figure  3-2  shows  a  three-dimensional  surface  plot  and  contour  plot  of  the 
gain  in  radiation  resistance  versus  ka  and  kR  for  six  air  bubbles,  at  1  atm  ambient 
pressure,  in  the  configuration  of  Figure  3-1.  The  gain  has  a  peak  value  of  1.95  at  ka  = 
0.0138  and  kR  =  2.87.  The  increasing  value  of  gain  as  kR  gets  small  is  not  interesting 
because  of  the  expected  divergent  behavior  [Ref.  4].  Inspection  of  the  gain  plot  reveals 
that  the  peak  in  gain  is  very  narrow  in  the  ka-di  recti  on.  The  reason  is  the  very  sharp 
peak  in  the  tuning  curve  for  a  resonant  bubble.  The  Q  of  a  resonant  bubble  is  about 
\lkar  [Ref.  8]. 

The  gain  in  radiation  resistance  was  also  plotted  for  a  bubble-like  passive 
acoustic  resonator  whose  resonant  value  of  ka  is  0.1.  This  is  shown  in  Figure  3-3.  An 
air  bubble  at  50  atm  hydrostatic  pressure  will  have  kar=0.\.  It  is  also  expected  that  a 
more  depth-independent  resonator  such  as  a  thin-shell  could  be  constructed  with  kar=0. 1 . 
The  Q  for  this  case  is  about  10  (as  compared  to  about  70  for  an  air  bubble  at  1  atm 
ambient  pressure.)  Figure  3-3  shows  a  much  broader  range  of  ka  values  with  substantial 
(>  1.5)  gain.  The  implication  of  these  observations  is  that  using  air  bubbles  at  1  atm 
ambient  pressure  would  not  be  practical,  but  that  using  resonators  with  kar=0. 1  could 
result  in  a  practical  application. 

Table  3-1  shows  numerical  results  for  the  planar  configuration  for  two 
through  ten  resonators  at  1  atm  ambient  pressure.  The  calculations  were  performed  using 
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Figure  3-2    Gain  in  radiation  resistance  versus  ka  and  kR  for  the  planar  configuration 
with  six  air  bubbles  at  1  atm  ambient  pressure  (a)  3D  surface  plot  (b)  contour  plot 
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Figure  3-3  Gain  in  radiation  resistance  versus  ka  and  kR  for  the  planar  configuration  for 
six  resonators  with  kar  =  0.1  (a)  3D  surface  plot  (b)  contour  plot 
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N 

ka 

kR 

Gain 

1  ird\ 

n  \ire,\ 

AIT1 

Dl 

2 

0.0138092 

2.65 

1.26 

0.386 

0.773 

90.0 

1.2 

3 

0.0138041 

2.69 

1.59 

0.573 

1.72 

90.0 

4.2 

4 

0.0138111 

2.82 

1.87 

0.636 

2.55 

90.0 

6.2 

5 

0.0138244 

2.86 

1.95 

0.553 

2.76 

90.0 

6.7 

6 

0.0138355 

2.87 

1.95 

0.465 

2.79 

90.0 

6.7 

7 

0.0138427 

2.87 

1.95 

0.398 

2.79 

90.0 

6.7 

8 

0.0138471 

2.87 

1.95 

0.349 

2.79 

90.0 

6.7 

9 

0.0138494 

2.87 

1.95 

0.310 

2.79 

90.0 

6.7 

10 

0.0138499 

2.87 

1.95 

0.279 

2.79 

90.0 

6.7 

Table  3-1    Planar  configuration  results  for  air  bubbles  at  1  atm  ambient  pressure. 


ka0 

ka 

kR 

Gain 

1  u*\ 

LIT1 

0.01 

0.0138355 

2.86783 

1.95490 

0.474090 

90.0 

0.025 

0.0138355 

2.86783 

1.95490 

0.473965 

90.0 

0.05 

0.0138355 

2.86785 

1.95488 

0.473517 

90.0 

0.075 

0.0138355 

2.86790 

1.95482 

0.472770 

90.0 

0.1 

0.0138355 

2.86801 

1.95472 

0.471723 

90.0 

0.2 

0.0138355 

2.86927 

1.95350 

0.464611 

90.0 

0.5 

0.0138342 

2.88818 

1.93501 

0.420116 

90.0 

0.75 

0.0138342 

2.92865 

1.89464 

0.367996 

90.0 

Table  3-2   Planar  configuration  results  for  six  resonators  allowing  ka0  to  vary. 
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Mathematica  with  a  precision  of  about  ten  significant  digits  (the  working  precision  was 
twenty  significant  digits).  Fewer  digits  are  shown  in  the  table  because  they  are  not  of 
general  interest.  Note  that  the  table  indicates  that  there  is  no  additional  increase  in  gain 
beyond  six  resonators.  This  is  in  part  the  basis  for  our  claim  that,  for  the  circular 
configurations,  six  resonators  is  the  best  choice.  Table  3-1  also  reveals  that  when  the 
gain  reaches  its  maximum  value,  the  total  volume  velocity  of  the  resonators  (#•  |  Uel  \ ) 
remains  constant  with  increasing  numbers  of  resonators.  The  phase  of  the  volume 
velocity  of  each  resonator  i/lF*)  is  90°  relative  to  the  transducer.  We  will  make  use 
of  this  fact  in  the  conical  configuration. 

Recall  that  the  calculations  for  Figures  3-1  and  3-2  and  Table  3-1  are  for 
ka0  =  0.2.    Table  3-2  shows  the  results  for  six  resonators  where  we  allow  ka0  to  vary. 


N 

ka 

kR 

Gain 

I  ir'\ 

N  \Urd\ 

LIT* 

Dl 

2 

0.0996243 

2.65 

1.26 

0.388 

0.776 

90.0 

1.2 

3 

0.0993435 

2.70 

1.59 

0.576 

1.73 

90.0 

4.2 

4 

0.0997272 

2.82 

1.87 

0.639 

2.56 

90.0 

6.2 

5 

0.100477 

2.86 

1.94 

0.556 

2.78 

90.0 

6.7 

6 

0.101117 

2.87 

1.95 

0.467 

2.80 

90.0 

6.7 

7 

0.101548 

2.87 

1.95 

0.400 

2.80 

90.0 

6.7 

8 

0.101813 

2,87 

1.95 

0.350 

2.80 

90.0 

6.7 

9 

0.101949 

2.87 

1.95 

0.311 

2.80 

90.0 

6.7 

10 

0.101981 

2.87 

1.95 

0.280 

2.80 

90.0 

6.7 

Table  3-3   Planar  configuration  results  for  resonators  with  kar=0A. 
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The  results  show  that,  over  a  large  range  of  values,  the  gain  in  radiation  resistance  is 
independent  of  the  value  of  ka$  chosen.  The  numerical  results  for  resonators  with 
kar=0. 1  are  shown  in  Table  3-3.  The  values  of  kR,  Gain,  IF1,  and  directivity  index  (D/) 
remain  the  same  as  in  Table  3-1 .  The  values  of  ka  are  different  but  it  can  be  shown  that 
the  ratio  kalka,  is  the  same.  This  indicates  that  over  a  large  range  of  values,  the  gain  in 
radiation  resistance  is  independent  of  the  value  of  kaT  for  bubble-like  resonators. 

3.       BEAM  PATTERNS  AND  DIRECTIVITY 

The  three-dimensional  beam  pattern  can  be  found  by  plotting  equation  2-24 
over  all  spherical  angles.    Figure  3-4  is  the  three-dimensional  far-field  beam  pattern  for 


Figure  3-4     Three-dimensional  beam  pattern  for  planar  configuration  with  six 
resonators  at  1  atm  ambient  pressure. 
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the  planar  configuration  with  six  air  bubbles  at  1  atm  ambient  pressure.  The  three- 
dimensional  (3D)  beam  pattern  provides  the  reader  with  an  idea  of  what  the  pressure 
field  will  look  like  in  3D  space  without  having  to  resort  to  multiple  2D  plots  in  various 
planes.  Now  that  we  know  what  the  beam  pattern  looks  like  we  can  better  utilize  our 
polar  plots  to  analyze  the  results.  From  the  3D  plot  alone,  we  can  see  that  there  is  some 
significant  directivity  along  the  acoustic  axis  (in  this  case  the  positive  z-axis).  We  can 
also  see  that  the  pressure  field  is  symmetric  about  the  z-axis  and  the  x-y  plane. 

The  usual  method  of  presenting  a  beam  pattern  is  the  dB  polar  plot.   Figure 
3-5  is  a  polar  plot  of  the  beam  pattern  for  the  planar  configuration  with  six  resonators. 


1 


A      \ 
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Figure  3-5    Polar  plot  of  the  beam  pattern  for  the  planar  configuration  with  six 
resonators  at  1  atm  ambient  pressure. 
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The  zero  dB  reference  is  the  pressure  that  would  exist  due  to  the  transducer  if  no  passive 
elements  were  present.  The  dashed  line  represents  the  sound  pressure  level  that  would 
exist  due  to  a  spherical  source  radiating  the  same  acoustic  power  as  the  transducer  does 
in  the  presence  of  the  passive  elements.  The  directivity  index,  defined  by  equation  2-26, 
can  be  read  directly  from  Figure  3-5  as  the  difference  in  dB  between  the  sound  pressure 
level  of  the  transducer  in  the  presence  of  passive  acoustic  radiators  and  that  of  the 
spherical  source  radiating  the  same  acoustic  power.  The  directivity  index  in  this  case  is 
6.7  dB.  Tables  3-1  and  3-3  list  the  directivity  index  for  the  planar  configuration  for  two 
through  ten  resonators  at  1  atm  ambient  pressure  and  for  kar=0A. 

B.      THE  CONICAL  CONFIGURATION 

1.       THE  GEOMETRY 

The  relative  volume  velocity  of  the  passive  acoustic  resonators  was  noted  in 
the  last  section  to  be  ninety  degrees  out  of  phase  with  the  transducer  volume  velocity. 
We  also  noted  that  the  far-field  beam  pattern  was  symmetric  about  the  x-y  plane.  For 
a  directional  source,  we  usually  want  the  sound  energy  to  be  directed  along  a  single 
acoustic  axis  with  any  side  or  back  lobes  minimized.  We  make  use  of  the  phase  of  the 
resonators'  relative  volume  velocity  in  the  conical  configuration  to  improve  the  directivity 
of  the  system.  In  the  conical  configuration  the  transducer  is  displaced  one-quarter 
wavelength  out  of  the  plane  of  the  resonators.  Figure  3-6  shows  the  geometry  of  the 
conical  configuration  for  the  case  of  six  passive  acoustic  resonators.  The  dimensionless 
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Figure  3-6   Geometry  of  the  conical  configuration  with  six  resonators. 

radius  of  the  circle  of  resonators  is  kR  and  the  dimensionless  height  of  the  transducer  is 

kH  =  tt/2. 

2.       GAIN  IN  RADIATION  RESISTANCE 

Surface  and  contour  plots  of  the  gain  in  radiation  resistance  for  the  conical 
configuration  with  six  resonators  at  1  atm  ambient  pressure  are  shown  in  Figure  3-7. 
The  plots  are  similar  to  the  planar  configuration  plots  except  that  the  gain  is  not  as  large. 
The  peak  in  gain  seen  in  Figure  3-7  is  1.82  at  &2=0.0138  and  kR=2.69.  The  value  of 
ka  is  1.005  times  the  resonant  value  for  a  single  bubble  at  1  atm  ambient  pressure  and 
is  approximately  the  same  as  the  value  of  ka  for  maximum  gain  in  the  planar 
configuration.    As  in  the  planar  configuration,  the  peak  in  gain  is  very  narrow   in   the 
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Figure  3-7   Gain  in  radiation  resistance  versus  ka  and  kR  for  the  conical  configuration 
with  six  resonators  at  1  atm  ambient  pressure  (a)  3D  surface  plot,  (b)  contour  plot 
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Figure  3-8    Gain  in  radiation  resistance  versus  ka  and  kR  for  the  conical  configuration 
with  six  resonators  for  kar  =  0.1  (a)  3D  surface  plot  (b)  contour  plot 


32 


/^-direction.  Figure  3-8  shows  the  surface  and  contour  plots  for  the  conical 
configuration  where  the  resonators  are  assumed  to  have  a  value  of  kar=0. 1 .  Once  again, 
the  peak  has  broader  support  in  the  ^-direction  for  kar=0A. 

The  numerical  results  for  the  conical  configuration  are  given  in  Tables  3-4, 
3-5,  and  3-6.   Table  3-4  gives  the  results  for  air  bubbles  at  1  atm  ambient  pressure  for 


TV 

ka 

kR 

Gain 

|  ud\ 

n  \irel  \ 

LUei 

Dl 

2 

0.0138079 

2.29 

1.29 

0.421 

0.841 

90.0 

4.3 

3 

0.0138038 

2.49 

1.58 

0.570 

1.71 

90.0 

6.7 

4 

0.0138115 

2.65 

1.78 

0.589 

2.36 

90.0 

8.1 

5 

0.0138238 

2.69 

1.82 

0.501 

2.51 

90.0 

8.4 

6 

0.0138330 

2.69 

1.82 

0.420 

2.52 

90.0 

8.4 

7 

0.0138384 

2.69 

1.82 

0.360 

2.52 

90.0 

8.4 

8 

0.0138409 

2.69 

1.82 

0.315 

2.52 

90.0 

8.4 

9 

0.0138413 

2.69 

1.82 

0.280 

2.52 

90.0 

8.4 

10 

0.0138399 

2.69 

1.82 

0.252 

2.52 

90.0 

8.4 

Table  3-4  Conical  configuration  results  for  air  bubbles  at  1  atm  ambient  pressure. 

two  through  ten  resonators.  The  conical  configuration,  like  the  planar  configuration, 
does  not  provide  any  additional  gain  in  radiation  resistance  for  more  than  six  resonators. 
Table  3-5  gives  the  results  for  six  resonators  in  a  conical  configuration  where  we  allow 
ka0  to  vary  from  the  usual  value  of  0.2.  The  gain  in  radiation  resistance  does  not  vary 
significantly  with  the  value  of  ka^  chosen.   Table  3-6  presents  the  results  for  the  conical 
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ka0 

ka 

kR 

Gain 

1  u*\ 

LU* 

0.01 

0.013833 

2.69149 

1.82404 

0.428152 

90.0 

0.025 

0.013833 

2.69149 

1.82402 

0.428040 

90.0 

0.05 

0.013833 

2.69151 

1.82404 

0.427641 

90.0 

0.075 

0.013833 

2.69155 

1.82403 

0.426980 

90.0 

0.1 

0.013833 

2.69165 

1.82402 

0.426062 

90.0 

0.2 

0.013833 

2.69275 

1.82393 

0.419965 

90.0 

0.5 

0.013833 

2.70910 

1.82189 

0.38413 

90.0 

0.75 

0.013833 

2.74319 

1.81422 

0.344514 

90.0 

Table  3-5   Conical  configuration  results  for  six  resonators  allowing  kaQ  to  vary. 


N 

ka 

kR 

Gain 

llf-l 

n  \  ir1  \ 

IIIrd 

DI 

2 

0.0995492 

2.29 

1.29 

0.423 

0.845 

90.0 

4.2 

3 

0.0993290 

2.49 

1.58 

0.573 

1.72 

90.0 

6.7 

4 

0.0997511 

2.65 

1.79 

0.592 

2.37 

90.0 

8.1 

5 

0.100443 

2.69 

1.82 

0.504 

2.52 

90.0 

8.4 

6 

0.100973 

2.69 

1.82 

0.422 

2.53 

90.0 

8.4 

7 

0.101288 

2.69 

1.82 

0.362 

2.53 

90.0 

8.4 

8 

0.101438 

2.69 

1.82 

0.317 

2.53 

90.0 

8.4 

9 

0.101459 

2.69 

1.82 

0.281 

2.53 

90.0 

8.4 

10 

0.101376 

2.69 

1.82 

0.253 

2.53 

90.0 

8.4 

Table  3-6     Conical  configuration  results  for  passive  acoustic  resonators  with 
ka  =0.1. 
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configuration  with  resonators  having  a  value  of  kar=0. 1.  The  values  of  kR,  Gain,  lfel, 
and  Dl  remain  about  the  same  as  in  Table  3-4.  While  the  values  of  ka  differ,  the  ratio 
kalkar  is  the  same  for  both  the  1  atm  bubble  and  the  kar=QA  resonator. 

Note  that,  for  four  or  more  resonators,  the  planar  configuration  has  a  larger 
gain  in  radiation  resistance  than  the  conical  configuration  with  the  same  number  of 
resonators.  Calculations  show  that  for  four  or  more  resonators,  the  maximum  gain 
always  occurs  when  the  transducer  is  in  the  plane  of  the  resonators. 

3.       BEAM  PATTERNS  AND  DIRECTIVITY 

Figure  3-9  shows  the  three-dimensional  far-field  beam  pattern,  defined  by 
equation  2-24  and  plotted  over  all  spherical  angles.  The  improvement  over  the  planar 
configuration  is  obvious.  We  traded  a  small  fraction  of  the  gain  in  radiation  resistance 
for  a  large  increase  in  directivity.  In  addition  to  concentrating  sound  energy  along  an 
acoustic  axis,  the  conical  configuration  also  has  smaller  side  and  back  lobes  than  the 
planar  configuration.  Figure  3-10  is  the  polar  dB  plot  of  the  conical  configuration  with 
six  resonators  at  1  atm  ambient  pressure.  From  Figure  3-10  we  note  that  the  conical 
configuration  concentrates  3  dB  more  sound  pressure  along  the  acoustic  axis  than  the 
planar  configuration.   This  represents  a  doubling  of  the  acoustic  intensity. 

The  directivity  index  is  determined  from  Figure  3-10  as  a  described  in  the  last 
section.  The  directivity  index  for  six  resonators  is  8.4  dB.  Tables  3-4  and  3-6  list  the 
directivity  index  for  two  through  ten  resonators  at  1  atm  ambient  pressure  and  for 
kar=0. 1.  The  directivity  index,  like  the  gain,  does  not  increase  with  more  resonators 
past  its  value  for  six  resonators.     The  high  gain  in  radiation  resistance  and  high 
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Figure  3-9  Three-dimensional  beam  pattern  for  the  conical  configuration  with  six 
resonators  at  1  atm  ambient  pressure. 

directivity  index  make  the  conical  configuration  the  best  choice  for  a  system  design.  As 
we  will  see  in  the  following  chapters,  the  3D  beam  pattern  for  the  conical  configuration 
is  the  most  desirable  of  the  configurations  considered. 


36 


Figure  3-10     Polar  plot  of  the  beam  pattern  on  a  dB  scale  for  the  conical 
configuration  with  six  resonators  at  1  atm  ambient  pressure. 
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IV.    LINE  CONFIGURATION  TYPE  I 

In  this  chapter  and  the  next  two,  we  extend  the  work  done  by  Ellsworth  by 
considering  linear  configurations  of  passive  acoustic  resonators.  The  motivation  is  to 
investigate  whether  a  linear  configuration  offers  any  improvement  in  power  radiated  and 
directivity  compared  to  a  circular  configuration.  We  present  results  for  six  resonators 
for  easy  comparison  to  the  results  in  the  last  chapter.  The  first  line  configuration 
considered,  termed  line  configuration  type  I,  consists  of  a  line  of  equally  spaced 
resonators  on  only  one  side  of  the  transducer.  In  the  next  chapter,  we  consider  an  equal 
number  of  resonators  on  either  side  of  a  transducer  in  a  line  (line  type  II).  In  Chapter 
VI,  we  will  consider  a  Yagi-Uda  [Ref.  5]  array  of  resonators  (line  type  III). 

A.      THE  GEOMETRY 

Figure  4-1  shows  the  line  configuration  type  I  with  six  resonators.  The  transducer 
is  at  the  origin  of  the  coordinate  system  and  the  resonators  are  along  the  positive  x-axis. 
The  dimensionless  distance  from  the  transducer  to  the  first  resonator  is  kL.  The 
dimensionless  distance  between  resonators  is  kl.  The  transducer  radius  is  ka0=0.2.  Each 
passive  acoustic  resonator  has  radius  kan,  where  the  subscript  n  indicates  the  individual 
resonator. 
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Figure  4-1    Geometry  of  the  line  configuration  type  I  with  six  resonators. 
B.      GAIN  IN  RADIATION  RESISTANCE 

Figure  4-2  shows  the  gain  in  radiation  resistance  presented  to  the  transducer  in  the 

presence  of  a  system  of  six  passive  acoustic  radiators  arranged  in  line  configuration  type 

I.    The  gain  in  radiation  resistance  is  plotted  versus  kL  and  kl  with  ka  held  constant  at 

the  value  that  provides  maximum  gain  for  a  finite  value  of  kL.  The  maximum  considered 

is  the  one  at  kL  =  3.6  and  kl=2.5,  for  which  ka=0. 01373.  The  gain  at  this  point  is  1.22 

(compared  to  the  conical  configuration  with  a  gain  of  1 .82  for  six  resonators).  The  peaks 

that  occur  in  the  gain  as  kL  goes  to  zero  were  only  briefly  considered.   This  is  discussed 

in  the  next  section.    A  plot  of  the  gain  in  radiation  resistance  versus  kL  and  kl  for 

ka=kar=0. 1  is  indistinguishable  from  Figure  4-2  and  is  not  presented  here.   We  expect 

that,  as  for  the  circular  configurations,  the  kar=0A  case  is  preferable  because  the  range 
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Figure  4-2   Gain  in  radiation  resistance  versus  kL  and  kl  for  the  line  configuration 
type  I  with  six  resonators  at  1  atm  ambient  pressure  (a)  3D  surface  plot  (b)  contour  plot 
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of  acceptable  values  of  ka,  i.e.  which  yield  a  gain  greater  than  a  certain  value,  would 
be  greater. 

Tables  4-1  and  4-2  present  results  for  the  line  configuration  type  I  for  one  through 
fourteen  resonators  at  1  atm  ambient  pressure  and  with  kar=0. 1  respectively.  Notice  that 
the  gain  remains  fairly  constant  at  about  1.2.  The  relative  volume  velocities  listed  are 
those  for  the  resonator  closest  to  the  transducer.  As  before,  the  primary  difference 
between  the  two  tables  is  the  values  of  ka.  It  can  be  shown  that  the  ratio  kalkar  is  the 
same  for  both. 

C.      BEAM  PATTERNS  AND  DIRECTIVITY 

Figure  4-3  shows  the  three-dimensional  far-field  beam  pattern  for  the  line  type  I 
configuration  with  six  resonators  at  1  atm  ambient  pressure.  The  acoustic  axis  is  along 
the  negative  x-axis.  Figure  4-4  shows  a  polar  plot  of  the  beam  pattern.  The  directivity 
index,  found  from  Figure  4-4  or  Table  4-1,  is  4.6.  The  line  configuration  type  I  can  be 
used  to  increase  the  power  output  and  directivity  of  an  underwater  transducer;  however, 
the  conical  configuration  results  are  clearly  better. 

In  the  last  section,  we  described  divergent  peaks  in  the  plot  of  gain  in  radiation 
resistance  for  vanishing  kL.  For  one  of  these  peaks,  for  example,  the  gain  2XkL  =  \  and 
&/=2.33  is  1.82.  We  plotted  a  beam  pattern  using  these  values  in  order  to  better 
understand  these  peaks  in  the  gain  plot.  The  directivity  index  was  found  to  be  DI=  1 1 .0. 
Figure  4-5  is  a  polar  plot  of  the  resulting  beam  pattern.  The  acoustic  axis  is  now  the 
positive  x-axis.   This  is  an  indication  that  the  peaks  as  kL  goes  to  zero  correspond  to  a 
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N 

*o 

kl 

kL 

Gain 

\TT*\ 

Z   IT" 

Dl 

1 

0.0138101 

— 

2.80 

1.11 

0.330 

90.0 

0.8 

2 

0.0137779 

2.09 

3.08 

1.16 

0.437 

86.2 

1.9 

4 

0.0137538 

2.32 

3.34 

1.20 

0.511 

86.1 

3.7 

6 

0.0137347 

2.49 

3.59 

1.22 

0.503 

86.0 

4.6 

8 

0.0137216 

2.62 

3.68 

1.22 

0.472 

87.3 

5.2 

10 

0.0137276 

2.60 

3.56 

1.24 

0.518 

86.6 

5.4 

12 

0.0137201 

2.64 

3.59 

1.24 

0.510 

90.5 

5.5 

14 

0.0138101 

2.67 

3.61 

1.25 

0.505 

91.1 

5.5 

Table  4-1     Line  configuration  type  I  results  for  air  bubbles  at  1  atm  ambient 
pressure. 


N 

ka 

kl 

kL 

Gain 

1  ^"1 

z  irl 

Dl 

1 

0.102413 

2.80 

1.11 

0.331 

90.0 

0.8 

2 

0.100606 

2.09 

3.08 

1.16 

0.439 

86.2 

2.0 

4 

0.0993406 

2.32 

3.34 

1.20 

0.514 

86.1 

3.7 

6 

0.0987282 

2.45 

3.46 

1.22 

0.530 

87.3 

4.6 

8 

0.0983271 

2.53 

3.52 

1.23 

0.528 

88.5 

5.1 

10 

0.0980380 

2.60 

3.57 

1.24 

0.522 

89.6 

5.4 

12 

0.0978263 

2.64 

3.59 

1.24 

0.514 

90.5 

5.5 

14 

0.0978263 

2.67 

3.61 

1.25 

0.509 

91.0 

5.5 

Table  4-2  Line  configuration  type  I  results  for  resonators  with  kar=0A. 
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Figure  4-3   Three-dimensional  beam  pattern  for  the  line  configuration  type  I  with 
six  resonators  at  1  atm  ambient  pressure. 

condition  on  the  values  of  kl  that  provide  constructive  interference.  It  is  probably  related 

to  the  total  length  of  the  aperture.    Further  investigation  is  required  to  fully  understand 

the  occurence  of  these  gain  peaks  and  to  determine  whether  a  practical  system  can  utilize 

this  phenomenon.  If  exploitable,  this  system  provides  the  best  directivity  of  any  analyzed 

in  this  research. 
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Figure  4-4    Polar  plot  of  the  beam  pattern  for  the  line  configuration  type  I  with  six 
resonators  at  1  atra  ambient  pressure 


...5 

V  '  '            o 

dB        """-W    '"■•. 

/A 

/     s  ' 

-  J X\) , 

/\    /NS,  ...-■»■ 

-  ^ 

v  / 

V               'J 

/  J  /     ■•' 

f         7^ V 

'*■-...          

.■'- LmeTyr*! 

SphcnaJ 

Figure  4-5    Polar  plot  of  the  beam  pattern  for  the  line  configuration  type  I  with  six 
resonators  at  1  atm  ambient  pressure  for  the  divergent  peak.  kL=\,  fc/=2.33 
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V.    LINE  CONFIGURATION  TYPE  II 

A.      THE  GEOMETRY 

Figure  5-1  shows  the  geometry  for  the  line  configuration  type  II  with  six 
resonators.  In  this  configuration,  half  of  the  resonators  are  on  one  side  of  the  transducer 
and  half  are  on  the  other  side.  The  dimensionless  distance  from  the  transducer  to  thefirst 
resonator  on  either  side  is  kL.  The  dimensionless  inter-resonator  spacing  is  kl.  The 
transducer  is  assumed  to  have  a  dimensionless  radius  of  ka0=Q.2. 


Jt 


Jl 


L 


A 


Figure  5-1    Geometry  of  the  line  configuration  type  II  with  six  passive  acoustic 
resonators. 
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B.      GAIN  IN  RADIATION  RESISTANCE 

The  gain  in  radiation  resistance  versus  kL  and  kl  for  a  fixed  ka  is  shown  in  Figure 
5-2.  The  value  chosen  for  ka  is  the  value  that  provides  maximum  gain  at  the  local 
maxima  at  kL=3.0  and  kl=2.2,  for  which  &2=0. 01376.  The  gain  at  this  point  is  1.56. 
As  with  the  line  type  I,  the  gain  plot  for  line  type  II  diverges  as  kL  goes  to  zero.  We 
have  not  investigated  this  behavior  in  detail.  Tables  5-1  and  5-2  give  the  numerical 
results  for  line  type  II  for  two  through  twenty  air  bubbles  at  1  atm  ambient  pressure  and 
for  resonators  with  kar=0A  respectively. 


N 

ka 

*/ 

kL 

Gain 

1  vd\ 

L  IT'1 

Dl 

2 

0.0138092 

2.65 

1.26 

0.386 

90.0 

0.7 

4 

0.0137761 

2.06 

2.87 

1.42 

0.631 

84.3 

2.7 

6 

0.0137595 

2.24 

2.99 

1.56 

0.752 

84.9 

2.8 

8 

0.0137478 

2.38 

3.06 

1.68 

0.852 

85.0 

4.9 

10 

0.0137387 

2.48 

3.10 

1.78 

0.923 

85.4 

5.5 

12 

0.0137315 

2.54 

3.17 

1.88 

0.979 

85.1 

5.9 

14 

0.0137247 

2.62 

3.15 

1.96 

1.014 

86.4 

6.1 

16 

0.0137192 

2.67 

3.16 

2.04 

1.046 

86.9 

6.3 

18 

0.0137143 

2.70 

3.17 

2.10 

1.073 

87.3 

6.4 

20 

0.0137099 

2.74 

3.18 

2.17 

1.096 

87.7 

6.5 

Table  5-1   Line  configuration  type  II  results  for  air  bubbles  at  1  atm  ambient 
pressure. 
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Figure  5-2   Gain  in  radiation  resistance  versus  Id  and  kl  for  the  line  configuration 
type  II  with  six  resonators  at  1  atm  ambient  pressure  (a)  3D  surface  plot  (b)  contour  plot 
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As  seen  in  Table  5-1,  the  gain  in  radiation  resistance  does  not  reach  a  constant  value 
(over  the  range  of  N  considered).  For  twenty  resonators  the  gain  is  2.17.  This  is  the 
highest  gain  calculated  in  this  research.  The  values  of  kl  and  kL  are  the  same  in  Tables 
5-1  and  5-2.  The  difference  is  in  the  values  of  ka.  It  can  be  shown  that  the  ratio  kalkar 
is  the  same  in  both  cases. 


N 

ka 

kl 

kL 

Gain 

1  ird\ 

L  Uel 

DI 

2 

0.102360 

2.65 

1.26 

0.388 

90.0 

0.7 

4 

0.100526 

1.89 

2.98 

1.42 

0.633 

84.2 

2.9 

6 

0.0996309 

2.24 

2.99 

1.56 

0.756 

84.9 

4.1 

8 

0.0990346 

2.38 

3.06 

1.68 

0.857 

85.1 

5.0 

10 

0.0985781 

2.48 

3.10 

1.78 

0.927 

85.5 

5.5 

12 

0.0982567 

2.53 

3.20 

1.87 

0.992 

84.7 

5.9 

14 

0.0978998 

2.62 

3.15 

1.96 

1.02 

86.4 

6.1 

16 

0.0976349 

2.66 

3.16 

2.04 

1.05 

86.9 

6.3 

18 

0.0974038 

2.70 

3.17 

2.11 

1.08 

87.3 

6.4 

20 

0.0971994 

2.74 

3.18 

2.17 

1.10 

87.6 

6.5 

Table  5-2   Line  configuration  type  II  results  for  resonator  with  kar=QA. 

C.      BEAM  PATTERNS  AND  DIRECTIVITY 

Figure  5-3  shows  the  three-dimensional  far-field  beam  pattern  for  the  line 
configuration  type  II  with  six  resonators  at  1  atm  ambient  pressure.  The  beam  pattern 
for  this  configuration  is  not  particularily  interesting.    The  acoustic  axis  is  along  the 
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Figure  5-3  Three-dimensional  beam  pattern  for  the  line  configuration  type  II  with 
six  resonators  at  1  atm  ambient  pressure. 

positive  and  negative  axes  and  there  is  some  directivity.  Figure  5-4  shows  the  dB  beam 
pattern.  The  directivity  index  (DF),  from  Figure  5-4  or  Table  5-1,  is  4.9  dB.  Unlike  the 
circular  configurations,  the  Dl  does  not  reach  a  constant  valve  (over  the  range  of  N 
considered).  For  twenty  resonators,  the  line  type  II  provides  a  directivity  index  of  6.5 
dB. 

Figure  5-5  shows  the  beam  pattern  for  this  configuration  when  using  the  value  of  kl 
corresponding  to  the  peak  at  kL—\.  The  value  of  kl  is  1.72,  the  directivity  index  is  8.0 
dB,  and  the  gain  is  2.15.     The  large  gain  and  high  directivity  makes  this  system 
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Figure  5-4    Polar  plot  of  the  beam  pattern  for  the  line  configuration  type  II  with  six 
resonators  at  1  atm  ambient  pressure 
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Figure  5-5    Polar  plot  of  the  beam  pattern  for  the  line  configuration  type  II  with  six 
resonators  at  1  atm  ambient  pressure  for  the  divergent  peak.  kL=\,  kl=U2 
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potentially  useful  for  practical  application.  However,  further  investigation  is  required 
to  understand  this  behavior  and  to  determine  whether  a  practical  system  can  utilize  the 
large  gain  in  radiation  resistance  at  small  values  of  kL.  We  provide  the  result  as  a  basis 
for  this  investigation. 
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VI.    AN  ACOUSTIC  YAGI-UDA  ARRAY 

The  Yagi-Uda  electomagnetic  antenna  system  [Ref.  5]  uses  parasitic  antenna 
elements  around  a  feed  antenna  element  in  order  to  direct  a  travelling  wave.  It  depends 
on  mutual  coupling  between  the  elements.  The  analogy  to  the  underwater  transducer 
with  passive  acoustic  resonators  is  direct.  The  transducer  is  the  feed  element  and  the 
passive  acoustic  resonators  are  the  parasitic  elements.  In  the  underwater  case,  the 
elements  are  coupled  by  the  sound  field.  We  are  basically  trying  to  do  the  same  thing 
that  antenna  farm  designers  have  been  doing  for  years  with  radio  waves.  Once  you  get 
past  the  obvious  analogy,  the  similarities  end.  The  acoustic  Yagi-Uda  array  described 
here  is  more  closely  related  to  the  work  described  earlier  in  this  paper  than  to  the 
electromagnetic  antenna  system.  We  will  alternately  refer  to  the  Yagi-Uda  array  as  line 
configuration  type  III.  In  the  sections  that  follow,  we  will  develop  the  theory  and  then 
present  numerical  and  graphical  results. 

A.      THEORY 

In  the  design  of  the  Yagi-Uda  configuration  (line  type  III),  we  are  interested  in 
directing  the  sound  energy  by  choosing  resonators  that  will  produce  constructive 
interference  in  one  direction  and  destructive  interference  in  the  opposite  direction.  We 
also  impose  the  condition  that  the  gain  in  radiation  resistance  be  not  less  than  one. 
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For  constructive  interference,  we  require  that  the  phase  of  the  pressure  wave 
outgoing  from  a  resonator  equal  the  phase  of  the  pressure  wave  incident  at  the  resonator 
from  the  transducer.  This  constraint  can  be  met  only  if  the  resonator  is  stiffness 
controlled.  We  then  place  an  additional  constraint  on  the  system  requiring  that  the  gain 
be  one  or  greater.  For  this  condition,  we  must  have  the  phase  of  the  pressure  wave 
incident  at  the  transducer  from  a  resonator  equal  to  the  phase  of  the  pressure  wave 
outgoing  from  the  transducer.  These  two  constraints  require  (1)  the  resonators  must  be 
stiffness-controlled,  and  (2)  the  dimensionless  distance  from  the  transducer  must  satisfy 

kdm    =    mir  ,  m  =  l,2,3,...  (6-1) 

where  d      =       the  distance  from  each  resonator  to  the  transducer 

and  k    =        wavenumber  in  the  medium. 

For  destructive  interference,  we  require  that  the  phase  of  the  pressure  wave 
outgoing  from  a  resonator  equal  x  plus  the  phase  of  the  pressure  wave  incident  at  the 
resonator  from  the  transducer.  This  constraint  can  only  be  met  if  the  resonator  is  mass- 
controlled.  We  then  use  the  constraint  that  the  gain  be  one  or  greater.  These  two 
contraints  require  (1)  the  resonators  must  be  mass-controlled,  and. (2)  the  dimensionless 
distance  from  the  transducer  must  satisfy 

kdm    =    -  +  mi    ,  w=l,2,3,...  (6-2) 

We  force  the  resonators  to  be  approximately  stiffness-controlled  yet  still  close  to 
resonance  by  setting  ka  =kar  -  \l2(kaf.    We  force  the  resonators  to  be  approximately 
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mass-controlled  and  likewise  close  to  resonance  by  setting  ka  =  kar  +  M2(karf.  These 
values  result  in  the  resonators  being  driven  at  a  frequency  which  is  one  resonant  half- 
width  above  and  below  their  resonant  frequency  (the  resonant  half- width  is  \ll(kar)2). 
Figure  6-1  is  a  diagram  of  the  line  type  III  configuration.  We  show  the  distances 
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Figure  6-1  Diagram  of  the  Yagi-Uda  configuration  (line  type  III)  for  six  resonators. 
in  terms  of  wavelengths  and  indicate  which  resonators  are  mass-  and  stiffness-  controlled. 
In  the  next  section  we  present  the  results  of  this  research. 

B.      RESULTS 

Computer  simulations  of  the  line  type  III  configuration  indicate  that  the  analysis  of 
the  theory  in  the  previous  section  is  correct.  We  did  not  perform  any  optimization  of 
parameters  on  this  configuration  as  we  did  for  previous  configurations.  The  results 
presented  here  are  for  the  distances  and  radii  described  in  the  previous  section. 

Figure  6-2  shows  a  three-dimensional  far-field  beam  pattern  for  a  line  type  in 
configuration  of  six  resonators.  The  acoustic  axis  is  along  the  positive  x-axis,  as 
expected.   Figure  6-3  shows  the  polar  plot  of  the  beam  pattern.   Table  6-1  lists  the  gain 
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N 

Gain 

Dl 

2 

0.83 

4.9 

4 

0.90 

5.9 

6 

0.94 

6.4 

8 

0.97 

6.5 

10 

0.99 

6.5 

Table  6-1  Results  for  the  Yagi-Uda  configuration  (line  type  III) 
in  radiation  resistance  and  directivity  index  for  two  through  ten  resonators.  The  gain  for 
six  or  more  resonators  is  about  one  and  the  directivity  index  is  about  6.5  dB.  Figure  6-4 
is  the  beam  pattern  and  Figure  6-5  the  polar  plot  for  the  case  of  ten  passive  acoustic 
radiators,  and  is  provided  for  comparison.  As  indicated  in  Table  6-1,  there  is  no  major 
difference  in  either  the  gain  in  radiaion  resistance  or  the  directivity  index  for  the  N  =  6 
and  N  =  10  cases. 
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Figure  6-2    Three-dimensional  beam  pattern  for  line  type  III  configuration  with  six  resonators 
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Figure  6-3  Polar  plot  of  the  beam  pattern  for  the  line  type  HI  configuration  with  six  resonators 
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Figure  6-4   Three-dimensional  beam  pattern  for  line  type  III  configuration  with  ten  resonators 
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Figure  6-5  Polar  plot  of  the  beam  pattern  for  the  line  type  HI  configuration  with  ten  resonators 
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vn.  DISCUSSION 

We  have  analyzed  five  configurations  of  resonators  and  provided  many  graphs  and 
tables  of  results.  Table  7-1  is  a  comparison  of  the  results  for  all  configurations.  The 
measures  of  effectiveness  that  we  chose  for  comparison  are  (1)  gain  in  radiation 
resistance,  (2)  directivity  index,  and  (3)  the  far-field  beam  patterns.   Table  7-1  lists  the 


N 

Planar 

Conical 

Line  Type  I 

Line  Type  II 

Line  Type  III 

Gain 

Dl 

dB 

Gain 

Dl 

dB 

Gain 

Dl 

dB 

Gain 

Dl 

dB 

Gain 

Dl 
dB 

2 

1.26 

1.2 

1.29 

4.3 

1.16 

1.9 

1.26 

0.7 

0.83 

4.9 

4 

1.87 

6.2 

1.78 

8.1 

1.20 

3.7 

1.42 

2.7 

0.90 

5.9 

6 

1.95 

6.7 

1.82 

8.4 

1.22 

4.6 

1.56 

2.8 

0.94 

6.4 

8 

1.95 

6.7 

1.82 

8.4 

1.22 

5.2 

1.68 

4.9 

0.97 

6.5 

10 

1.96 

6.7 

1.82 

8.4 

1.24 

5.4 

1.78 

5.5 

0.99 

6.5 

Table  7-1    Comparison  of  the  results  for  all  configurations 

data  for  the  first  two  measures  of  effectiveness  and  Figure  7-1  provides  a  comparison  of 
the  far-field  beam  patterns  on  a  single  plot.  Recall  that  0  dB  represents  the  transducer 
with  no  resonators  present.  We  considered  the  case  of  six  resonators  throughout  the 
investigation  as  a  basis  for  comparison  of  the  three-dimensional  beam  patterns.  The 
N  =  6  case  is  highlighted  in  Table  7-1.  The  planar  configuration  has  the  best  gain,  almost 
double,  and  the  conical  configuration  has  the  best  directivity  index,  8.4  dB.   Figure  7-1 
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Figure  7-1    Polar  plot  of  the  far-field  beam  pattern  for  all  configurations  for  six 
resonators. 
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shows  that  the  conical  configuration  also  has  the  best  front-to-back  discrimination  along 
the  acoustic  axis.  Of  the  results  presented,  clearly  the  conical  configuration  is  the  best 
choice  for  a  real  system.  The  linear  configurations  are  generally  worse  than  the  circular 
configurations,  but  the  line  type  III  configuration  is  of  particular  interest.  The  addition 
of  six  passive  acoustic  resonators  in  a  Yagi-Uda  antenna-like  configuration  greatly 
improves  the  directivity  of  the  transducer.  However,  since  the  passive  resonators  are  not 
driven  at  their  resonant  frequency,  this  configuration  does  not  provide  any  power  gain. 

Further  investigation  of  the  divergent  peaks  in  gain  for  the  line  configurations  type 
I  and  II  is  warranted.  The  gain  diverges  as  the  distance,  L,  to  the  first  resonator  tends 
toward  zero.  It  may  be  possible  to  exploit  regions  of  finite  kL  for  which  this  behavior 
occurs.  Certain  values  of  resonator  spacing  satisfy  a  resonant  condition  related  to  the 
total  length  of  the  aperture  which  apparently  results  in  substantial  gain.  We  presented 
beam  patterns  for  a  system  operating  on  one  of  these  divergent  gain  peaks  and  found  that 
the  directivity  and  gain  of  the  system  were  better  than  any  other  system  analyzed. 
Further  investigation  into  this  behavior  is  encouraged. 

We  have  shown  that  resonators  with  a  value  of  kar  =  0. 1  provide  identical  results 
in  terms  of  gain  and  directivity  as  do  bubbles  at  1  atm  ambient  pressure.  We  also 
demonstrated  that  these  resonators  are  not  as  sensitive  to  the  value  of  ka,  that  is, 
substantial  gain  (>1.5)  may  be  achieved  over  a  reasonable  range  of  operating 
frequencies.  This  is  understood  as  follows.  For  an  air-filled  bubble  at  1  atm  ambient 
pressure,  Q  =  \/kar  is  about  70.  For  a  system  to  be  feasible,  the  dimensions  of  each 
resonator  must  be  matched  to  substantially  better  than  one  part  in  Q.   It  is  probably  not 
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reasonable  to  expect  that  a  collection  of  bubbles  can  be  matched  closely  enough  in  radius 
for  a  system  at  1  atm  ambient  pressure  to  be  practical.  It  may  be  possible  to  build, 
however,  a  system  of  bubble-like  resonators  such  as  thin-walled  shells  that  will  have  a 
kar  of  about  0. 1  and  therefore  a  Q  of  about  10. 

Table  7-2  describes  the  physical  size  of  a  system  for  various  configurations  of  six 
resonators  as  described  in  this  thesis.  The  overall  dimensions  of  the  system  are  the  same 


Frequency: 

10  Hz 

100  Hz 

1  kHz 

Radius  of  (cm): 
Resonator 

33 

240 

3.3 

24 

0.33 

2.4 

Aperture  size  of  (m): 
Planar 

137 

13.7 

1.37 

Conical 

128 

12.8 

1.28 

Line  Type  I 

383 

38.3 

3.83 

Line  Type  II 

357 

35.7 

3.57 

Line  Type  III 

413 

41.3 

4.13 

Table  7-2    Physical  size  of  the  system  of  resonators  for  10  Hz,  100  Hz,  and  1  kHz  for 
various  configurations  of  six  resonators. 

for  bubbles  at  1  atm  ambient  pressure  and  for  Jcar=0.l  resonators.  The  size  of  the 
resonator  increases  for  larger  values  of  kar.  The  conical  configuration,  chosen  as  the 
best  in  terms  gain  and  directivity,  has  the  smallest  physical  dimensions.  The  aperture 
size  (diameter  for  circular  and  total  length  for  linear  configurations)  for  the  conical 
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configuration  at  100Hz  is  about  13  meters.    A  kar  =  OA   resonator  for  this  case  would 
have  a  radius  of  24  cm. 
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vm.   CONCLUSION 

A.  SUMMARY 

We  developed  the  equations  for  the  network  analysis  of  a  transducer  in  the  presence 
of  passive  acoustic  radiators.  These  equations  were  implemented  in  computer  programs 
and  applied  to  five  configurations  of  resonators.  We  first  considered  the  two 
configurations  previously  examined  by  Ellsworth.  We  provided  results  for  the  case  of 
six  resonators  and  developed  three-dimensional  beam  patterns.  We  then  examined  two 
linear  configurations  and  provided  results.  Next,  we  developed  the  idea  of  a  Yagi-Uda 
antenna-like  array  of  passive  acoustic  resonators.  Finally,  we  conducted  a  comparison 
of  all  five  configurations. 

B.  CONCLUSIONS 

We  found  that  the  conical  configuration  of  six  resonators  is  the  best  choice  in  terms 
of  the  gain  in  radiated  power  output,  directivity  along  an  acoustic  axis,  and  minimum 
required  number  of  resonators.  A  linear  configuration  could  be  utilized  to  improve  the 
gain  and  directivity  of  a  transducer,  but  not  as  well  as  the  conical  configuration.  We 
found  that  the  acoustic  Yagi-Uda  antenna  configuration  can  be  used  to  improve  the 
directivity  of  an  underwater  transducer. 
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C.      RECOMMENDATIONS 

The  results  of  the  research  presented  in  this  paper  are  promising.  We  recommend 
that  further  work  be  done  to  (1)  design  a  resonator  with  kar  =  0.1,  (2)  build  a  system 
of  resonators  and  verify  the  theory  and  results  in  this  paper,  (3)  study  the  Yagi-Uda 
configuration  in  more  detail  including  optimizations  to  improve  the  gain  ,  and  (4)  study 
the  divergent  peaks  in  gain  as  kL  goes  to  zero  for  the  linear  configurations  to  determine 
whether  a  practical  system  can  be  designed  to  utilize  the  higher  theoretical  gain  and 
directivity. 
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APPENDIX 

The  following  Mathematica  program  code  is  an  example  of  how  programs  were 
written  to  perform  the  calculations  in  this  paper.  The  example  consists  of  the  standard 
definitions  for  the  conical  configuration  and  code  to  create  plots  as  found  in  this 
document  and  calculate  the  information  found  in  tables  for  the  conical  configuration. 
The  lines  in  bold  face  type  are  comments. 

Distance,  Radius,  and  Area 

ClearAll[ka,S,kl,krange] 

ka[0]  :=  kaO 

ka[nj  :=  kai 

S[nJ  :  =  ka[n]  A  2 

kl[n_,nj  :=  0 

kl[n_,0]  :=  Sqrt[kR"2  +  kH"2] 

kl[0,nj  :=  Sqrt[kRA2  +  kJT2] 

kl[n_,mj  :=  kR  Sqrt[2-2*Cos[(2  Pi/Num)  (n-m)]] 

krange[0,the_,ph_]  :  =  kr 

krange[n_,the_,phj  :=  Sqrt[ 

(kr  Sin[the]  Cos[ph]  -  kR  Cos[(2  Pi  n)/Num])A2  + 

(kr  Sin[the]  Sin[ph]  -  kR  Sin[(2  Pi  n)/Num])A2  + 

(kr  Cos[the]  -  kH)"2] 

Impedances 

ClearAll[Za,Zar,Dfc] 
Dfc[nJ  :=  Exp[j  ka[n]]  /  (l+ka[n]) 
Za[n  J  :  =  -  I  (S[0]/S[n])  ((3  gam  Po)/(p  cA2  ka[n])) 
Zar[n_,nJ  :=  (S[0]/S[n])  (I  ka[n]/(l+I  ka[n])) 
Zar[n_,m_]  :=  Dfc[n]  Zar[m,m]  (ka[m]/kl[n,m])  Exp[ 
I  (ka[m]  -  kl[n,m])] 
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Gain  in  radiation  resistance 

Clear  All  [ZB,ZC,  Gain] 
ZB[nJ  :=  Zar[n,0] 
ZC[n_,nJ  :=  Zar[n,n]  +  Za[n] 
ZC[n_,m_]  :=  Zar[n,m] 

Gain  :=  Re[Zar[0,0]  -  Array [ZB,Num]  .  Inverse[Array[ 
ZC,{Num,Num}]]  .  Array[ZB,Num]]/Re[Zar[0,0]] 

Relative  velocities  and  field  pressure 

Clear[theta,phi,the,ph,n] 

ClearAll[BP,dBBP,Prad,Urel] 

Urel  :=  -  Inverse[Array[ZC,{Num,Num}]]  .  Array [ZB,Num]; 

Prad[0,the_,phj  :=  Zar[0,0]  (ka[0]/krange[0,the,ph])  Exp[ 

I  (ka[0]-krange[0,the,ph])] 
Prad[n_,the_,ph_]  :=  Zar[n,n]  Urel[[n]]  (ka[n]/ 

krange[n,the,ph])  Exp[I  (ka[n]-krange[n,the,ph])] 
BP[the_,phJ  :=  Abs[Sum[Prad[n,the,ph],{n,0,Num}]]/ 

Abs[Prad[0,the,ph]] 
dBBP[the_,phJ  :=  20  Log[10,BP[the,ph]] 

Constants 

Clear[kH,kai,kaO,kR,kr,kar,gam,Po,c,p] 

gam  =  1.402; 

Po  =  1.013  10A5; 

c  =  1500; 

p  =  998; 

kar  =  Sqrt[(3  gam  Po)/(p  cA2)] 

Num  =  6 

kH=Pi/2 

ka0=0.2 

Create  surface  plot  of  Gain  versus  ka  and  kR 

plotl=Plot3D[Gain,{kR,l,5},{kai,0.9kar,l.l  kar}, 
PlotPoints-  >  40,PlotRange-  >  {0,2}, 
AxesLabel- >  {"kRVkaY'Gain"}] 
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Create  contour  plot 

plot2  =  Show[ContourGraphics[plotl], 
Frame- >  True, 
FrameLabel->  {"kR'V'ka"}, 
RotateLabel-  >  True, 
ContourSmoothing-  >  3, 
ContourShading-  >  False, 
AspectRatio- >  1, 
Contours- >  {0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.8}] 

Find  Maximum  Gain 

Clear  [kai,kR] 
{maxgain,{kaiRule,kRrule}}  = 

FindMinimum[-Gain,{kai,{kar,  1.002  kar}}, 
{kR,{2.7,2.8}}] 
kai  =  kai/.kaiRule 
kR=kR/.kRrule 
MaxGain  =  -maxgain 

Find  relative  volume  velocities 

N[{Abs[Urel[[l]]],Arg[Urel[[l]]]  180/Pi}] 

Find  directivity  index 

thAxis  =  0;  (*  theta  on  acoustic  axis  *) 

phAxis  =  0;  (*  phi  on  acoustic  axis     *) 

DI  =  dBBP[thAxis,phAxis]  - 

20  Log[10,Sqrt[MaxGain]]  //N 

Create  3D  beam  pattern 

BP2[th_,ph_]  =  Evaluate[BP[th,ph]]; 
Needs["Graphics'ParametricPlot3D'"] 
plot3=ParametricPlot3D[BP2  [theta,  phi] 

{Sin[theta]  Cos[phi], 

Sin[theta]  Sin[phi], 

Cos[theta]}, 

{theta,0,  Pi,  Pi/40}, 

phi,0,  2  Pi,  Pi/20}] 
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Create  2D  (dB)  beam  pattern  in  z-x  plane 

data=Table[N[dBBP[theta,0]], 

{theta,0,2  Pi, Pi/40}]; 
mystuff=Table[(data[[i]]  +  12) 

(N[Sin[(i-l)  Pi/40]], N[Cos[(i-l)  Pi/40]]}, 
{i,l,81}]; 
myplot=ListPlot[mystuff,PlotJoined->True, 
AspectRatio-  >  Automatic] 

The  following  requires  a  separately  created  template  and  label: 
Show[template, label,  myplot, 
AspectRatio-  >  Automatic, 
Axes-  >  False,  PlotRange-  >  All] 
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